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I .  Introduction 

This  report  contains  user  instructions,  a  listing  and 
documentation  for  a  microcomputer  BASIC  program  that  can  be  used 
to  compute  an  estimate  of  the  probability  that  an  electrometer 
based  detection  system  will  detect  a  submarine  electric  dipole 
field  during  an  encounter. 

The  program  generates  detection  probabilities  based  on  two 
encounter  models.  In  the  first  encounter  model,  the  detection 
system  uses  a  crosscorrelation  detector.  In  the  second  encounter 
model,  the  detection  system  uses  a  square  law  detector.  Relative 
to  operationally  realizable  values,  probabilities  based  on  the 
first  model  represent  upper  bounds  and  those  based  on  the  second 
model  represent  lower  bounds.  For  both  encounter  models,  the 
signal  is  proportional  to  the  square  of  the  component  of  the 
electric  field  that  is  measured  by  a  single  axis  electrometer  in 
an  infinite  medium  of  constant  conductivity  and  the  electric 
field  noise  does  not  change  with  changes  in  the  position  of  the 
electrometer.  Also,  for  both  encounter  models,  the  encounter  is 
a  straight  line  encounter  with  constant  vertical  separation. 

The  encounter  models  can  be  interpreted  as  models  of  a 
electric  anomaly  detection  system  that  is  moving  with  constant 
course,  speed  and  altitude  in  an  encounter  with  a  submarine 
moving  with  constant  course,  speed  and  depth.  Or,  they  can  be 
interpreted  as  models  of  a  stationary  electric  anomaly  detection 
system  in  an  encounter  with  a  submarine  moving  with  constant 
course,  speed  and  depth. 


The  program  parameters  include:  electrometer  axis 
direction,  electrometer  axis  depression  angle,  encounter  medium 
conductivity,  submarine  electric  moments,  submarine  course,  speed 
and  altitude,  electrometer  course,  speed  and  depth,  encounter 
lateral  range  (the  horizontal  range  at  the  closest  point  of 
approach  in  a  straight  line  encounter)  and  false  alarm  rate.  (In 
the  program,  a  false  alarm  is  the  event  that  the  detection  system 
classifies  noise  as  a  dipole  signal.) 
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II.  Program  User  Instructions 

As  listed  in  Appendix  10,  the  program  can  be  run  under  a 
BASIC  language  that  is  compatible  with  IBM  PC  BASIC.  If  the 
listing  is  used  to  enter  the  program  through  the  keyboard,  then 
the  program  should  be  saved  with  the  name  EAD.BAS  or  the  value  of 
N$  on  line  40  should  be  changed  to  the  file  name  under  which  it 
is  saved. 

The  program  contains  user  instructions  in  the  form  of  query 
and  parameter  limitation  messages.  As  an  example  of  the  former, 
after  starting  the  program  under  BASIC,  the  following  message 
should  appear: 

Electric  Anomaly  Detection  (EAD)  Lateral  Range  Function 

generate  or  print  a  program  data  file  (g/p)? 

By  entering  p,  data  can  be  printed  from  a  program  data  file 
that  was  generated  by  the  program.  By  entering  g,  a  program 
data  file  can  be  generated  for  a  set  of  user  specified 
conditions.  With  either  response,  a  sequence  of  additional 
queries  is  displayed.  These  queries  require  either  an  indicated 
response  or  a  parameter  value  as  the  input.  If  the  initial 
response  is  g,  the  sequence  includes  queries  whose  responses 
determine  whether  or  not  an  auxiliary  data  file  will  be  generated 
that  can  be  used  for  future  input  of  magnetic,  processing  or 
kinematic  data.  In  particular,  the  first  query  in  the  sequence 
gives  the  option  of  using  a  combined  electric,  processing  and 
kinematic  data  file.  If  the  response  indicates  that  it  should 
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be  used  and  the  file  is  available,  the  parameter  values  that 
remain  to  be  entered  in  order  to  generate  a  program  data  file  are 
the  following:  the  false  alarm  rate,  the  electric  field  noise, 
the  maximum  encounter  lateral  range  and  the  lateral  range  step. 
The  combined  file  should  be  used  only  if  the  effect  of  varying 
just  one  or  more  of  these  parameters  is  desired.  If  the  response 
indicates  that  the  file  will  not  be  used,  queries  concerning 
electric,  kinematic  and  processing  parameter  values  are 
displayed. 

After  all  of  the  program  parameter  values  have  been 
entered,  a  query  is  displayed  giving  the  option  of  generating  the 
combined  file.  Then  a  query  giving  the  option  of  generating  a 
program  data  file,  a  query  giving  the  option  of  printing  the 
encounter  parameter  values  and  a  query  giving  the  option  of 
printing  lateral  range  function  values  are  displayed.  The 
lateral  range  function  values  are  the  encounter  detection 
probabilities  indexed  by  lateral  range.  The  parameters  maximum 
lateral  range  and  lateral  range  step  determine  the  index  lateral 
ranges  of  the  encounters  for  which  probabilities  of  detection  are 
computed . 

The  program  generates  electric  signal  values  that  correspond 
to  points  in  time  during  an  encounter.  Following  the  lateral 
range  function  query,  a  query  is  displayed  that  gives  the  option 
of  printing  the  electric  signal  values  for  an  encounter.  If  the 
option  is  exercised,  the  option  is  repeated.  When  the  option  is 
not  exercised,  a  query  is  displayed  that  gives  the  option  of 
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generating  or  printing  a  new  program  data  file.  If  this  option 
is  not  exercised,  the  program  ends. 

Some  suggested  guides  for  determining  parameter  values  can 
be  found  in  Section  IV  of  this  report. 
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III.  Encounter  Model  Limitations 

In  the  two  encounter  models,  the  signal  is  the  unfiltered 
signal  that  is  generated  by  an  electric  dipole  that  moves 
relative  to  an  electrometer.  Depending  on  the  input  filtering, 
describing  a  submarine  electric  anomaly  field  as  an  electric 
dipole  field  should  not  be  a  significant  limitation  for  an 
encounter  slant  range  at  the  closest  point  of  approach  (CPA)  that 
is  greater  than  one  hull  length.  The  detection  decision  is  based 
on  samples  from  a  single  time  interval  (window)  that  is  centered 
on  the  CPA.  The  length  of  the  sampling  interval  and  the  sampling 
rate  are  parameter  values  that  are  inputs  to  the  program.  In 
terms  of  signal-to-noise  ratio,  there  is  an  optimum  sampling 
interval  length  (integration  time)  and  sample  rate.  Although 
dipole  signal  energy  is  not  symmetrically  distributed  about  the 
CPA  time,  for  a  given  sampling  interval,  the  difference  between 
the  signal  energy  for  the  optimum  interval  location  and  the  CPA 
centered  location  should  not  be  significant  in  most  cases. 

The  encounter  model  electrometer  noise  samples  are 
determined  by  a  gaussian  random  process  and  they  are  values  of 
independent  identically  distributed  normal  random  variables.  The 
standard  deviation  of  these  random  variables  is  referred  to  as 
the  electric  noise  and  the  variance  is  the  noise  in  the  sense  of 
the  signal-to-noise  ratio.  The  degree  of  correspondence  between 
this  process  and  operational  noise  depends  on  the  nature  of  the 
dominant  operational  noise  sources  and  on  the  electrometer  input 
fi^ering  (noise  whitening). 
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In  the  encounter  models,  the  intervals  are  adjacent  but  not 
overlapping  and  a  detection  statistic  corresponds  to  each  sample 
interval  and  its  value  is  determined  by  the  sample  values.  If 
the  value  of  the  detection  statistic  for  an  interval  equals  or 
exceeds  a  threshold  value,  a  detection  is  indicated.  The 
threshold  value  is  determined  by  the  false  alarm  probability 
which  in  turn  is  determined  by  the  false  alarm  rate  and  the 
sampling  interval  length. 

A  detector  that  used  a  moving  sample  interval  that  was 
generated  by  replacing  the  oldest  sample  by  the  newest  one  would 
correspond  more  closely  to  the  detector  in  an  operational 
detection  system.  In  a  model  of  the  detector,  since  the  sample 
windows  overlap,  the  detection  statistics  would  represent  a 
sequence  of  random  variables  that  were  correlated  over  an 
interval  equal  to  the  width  of  the  sample  interval.  Because  of 
this  dependence,  it  seems  unlikely  that  the  results  that  would  be 
obtained  with  an  encounter  model  based  on  an  overlapping  interval 
detector  would  differ  significantly  from  those  obtained  with  an 
encounter  model  based  on  a  nonoverlapping  interval  detector. 

For  encounter  lengths  of  the  order  of  a  few  nautical  miles 
or  less,  the  straight  line  encounter  condition  should  not  be  a 
significant  limitation.  In  particular,  this  should  be  the  case 
for  a  fixed  electrometer  since,  for  a  submarine  target,  vertical 
separation  and  course  changes  should  be  less  likely  to  occur. 

Other  models  are  available  that  can  be  used  as  the  basis  for 
computing  an  estimate  of  the  probability  that  an  electric  anomaly 
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detection  system  will  detect  a  submarine  during  an  encounter. 

For  example,  one  that  is  described  in  Appendix  6  could  be  used  to 
determine  the  slant  ranges  of  straight  line  encounters  for  which 
the  detection  probability  is  equal  to  a  specified  value.  The 
parameter  values  that  are  required  to  do  this  are  an  average 
submarine  electric  dipole  moment,  a  detection  system  capability 
factor  and  a  noise  factor.  Values  for  these  parameters  can  be 
determined  from  operational  data.  However,  the  values  are 
specific  to  averages  over  a  particular  set  of  encounter 
conditions.  An  advantage  of  the  two  encounter  models  relative  to 
this  model  is  their  adaptability  to  different  electric, 
processing  and  kinematic  conditions. 
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IV.  Parameter  Values 

The  electric  parameter  value  queries  are  generally 
explanatory  with  regard  to  the  value  that  should  be  entered. 

This  is  also  true  of  the  kinematic  parameter  value  queries. 
However,  there  is  some  ambiguity  with  respect  the  processing 
parameter  value  queries  and  the  noise  parameter  value  query.  To 
reduce  this  ambiguity,  a  brief  discussion  of  the  common 
characteristics  of  the  two  encounter  model  processing  parameters 
and  the  noise  parameter  is  given  below.  This  is  followed  by 
some  guidelines  for  choosing  these  parameter  values. 

In  both  encounter  models,  a  decision  is  made  at  the  end  of 
each  sampling  interval.  The  decision  is  either  noise  energy  was 
present  during  the  interval  or  noise  energy  plus  signal  energy 
was  present  during  the  interval.  The  sample  intervals  are 
adjacent,  equal  width,  nonoverlapping  time  intervals.  The  number 
of  samples  that  are  input  in  a  sample  interval  is  determined  by 
the  sampling  rate  and  the  interval  length. 

The  program  default  choice  for  the  sampling  rate  is  2  MAXF 
where  MAXF  is  a  parameter  that  is  labeled  the  maximum  electric 
signal  frequency.  This  sample  rate  is  the  Nyquist  rate  for  an 
ideal  low  pass  filter.  However,  the  signal  in  a  sample  interval 
that  is  computed  by  the  program  represents  an  unfiltered  dipole 
signal.  This  is  a  reasonable  approximation  if  the  signal  energy 
that  is  associated  with  signal  components  greater  than  MAXF  is 
relatively  small.  As  discussed  below,  the  noise  energy  in  a 
sample  interval  should  be  considered  to  be  proportional  to  MAXF 
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in  order  to  be  consistent  with  the  encounter  models.  Ideally,  a 
default  choice  for  MAXF  would  make  the  ratio  of  the  signal  energy 
to  the  noise  energy  a  maximum  for  the  sampling  interval  of  an 
encounter.  The  program  default  choice  for  MAXF  is  2 • MAXVM/MINRO 
where  MAXVM  is  a  user  estimated  maximum  encounter  relative  speed 
converted  from  knots  to  meters  per  second  and  MINRO  is  a  user 
estimated  maximum  slant  range  at  CPA  in  meters  in  terms  of  a  just 
detectable  target.  (A  precise  definition  of  a  just  detectable 
target  can  be  made  in  terms  of  a  specified  detection  probability, 
false  alarm  probability  and  target  dipole  moment.)  The  default 
choice  for  MAXF  is  consistent  with  the  observation  in  Reference  1 
that  if  an  optimum  value  for  MAXF  is  determined  for  a  minimum 
dipole  moment  target,  then  no  significant  increase  in  MAXF  is 
required  in  order  to  maintain  a  required  detection  probability  if 
the  encounter  lateral  range  is  decreased  even  though  the  signal 
energy  spectrum  is  shifted  to  higher  frequencies. 

Because  of  the  detection  statistics  and  the  gaussian  noise 
model,  if  there  were  no  penalty  for  decision  delay,  a  sample 
interval  length  for  a  signal  should  be  chosen  equal  to  the  signal 
duration,  since  this  would  make  the  detection  probability  for  an 
encounter  a  maximum.  In  the  program,  a  default  choice  for  the 
sample  interval  length  (integration  time)  is  2 • MAXRO/MINVM  where 
MINVM  is  a  user  estimated  maximum  encounter  relative  speed 
converted  from  knots  to  meters  per  second  and  MAXRO  is  a  user 
estimated  maximum  slant  range  at  CPA  in  meters  for  a  detectable 
encounter  (in  terms  of  a  specific  detection  probability  and  false 
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alarm  probability)  with  a  user  specified  maximum  dipole  moment 
target.  This  choice  might  be  considered  a  balance  between 
minimizing  decision  delay  and  maximizing  detection  probability. 

In  the  encounter  models,  the  sample  intervals  are  located  so  that 
the  CPA  time  is  at  the  center  of  a  sample  interval  and  this 
sample  interval  is  the  only  one  that  contains  signal  energy. 

This  characteristic  is  consistent  with  the  default  choice  for  the 
sample  interval.  For  a  given  sample  interval  length,  although  in 
general  the  interval  location  is  not  the  optimum  one  in  terms  of 
signal  energy,  it  should  be  approximately  so  in  most  cases. 

The  program  noise  parameter  is  SIG.  It  represents  the 
standard  deviation  a  associated  with  the  electric  noise  process 
of  the  two  encounter  models.  In  terms  of  the  ideal  low  pass 
filter  implied  by  the  encounter  models,  its  square  should  be 
equal  to  MAXF(SIGO)*  where  SIGO  is  the  electric  noise  process 
constant  power  spectral  density.  The  program  does  not  enforce 
this  relation.  Therefore,  in  using  the  program,  the  implied 
relation  between  the  two  input  parameters:  electric  noise  and 
maximum  electric  signal  frequency  should  be  kept  in  mind.  If  an 
average  value  of  the  peak-to-peak  electric  noise  for  an  encounter 
can  be  estimated,  for  example  from  a  electrometer  trace,  then  the 
value  for  SIG  should  be  chosen  so  that  the  estimate  is  4  to  6 


times  this  value. 


Appendix  1.  The  Detection  Statistics 

In  this  appendix,  y^,  y2,  ' ' ‘  /  Ym  are  sequential  values 
(voltages)  representing  the  sample  values  in  a  sample  interval. 
They  are  the  input  to  an  electrometer’s  detector.  With  these  m 
sample  values,  the  detector  computes  the  value  of  a  detection 
statistic.  This  value  is  represented  by  x  and  the  detection 
decision  corresponding  to  the  sample  interval  is  determined  by 
the  decision  rule:  If  x  >  x*,  then  the  input  during  the  sample 
interval  was  noise  plus  signal,  otherwise,  the  input  was  noise. 

For  both  encounter  models,  the  detection  probability  and  the 
false  alarm  probability  are  decreasing  functions  of  x*  and  the 
relation  is  one-to-one  in  both  cases.  In  the  program,  the  false 
alarm  probability  pf  is  used  to  determine  a  unique  value  of 
the  threshold  x*.  This  value  is  then  used  to  determine  a  unique 
value  of  the  detection  probability  p<j. 

In  the  program,  pf  is  found  using  the  relation  pf  =  R-6t 
where  R  is  the  false  alarm  rate  in  false  alarms  per  second 
and  6t  is  the  sample  interval  length  in  seconds.  This  relation 
is  based  on  the  following  argument:  With  no  signal  energy  in  a 
sample  interval,  yx  =  nlf  y2  =  n2,  /  ym  =  nm  where  n1#  n2, 

,  nm  are  noise  values  (voltages)  input  to  an  electrometer's 
detector.  The  n^  are  values  of  independent  normal  (gaussian) 
random  variables,  each  with  mean  zero  and  standard  deviation  a. 
Because  of  this,  in  the  encounter  models,  values  of  x  for 
different  sample  intervals  are  the  values  of  independent  random 
variables  that  determine  two  outcomes:  x  >  x*  or  x  <  x*. 


Therefore,  in  terms  of  these  outcomes,  an  encounter  is  a  sequence 
of  independent  Bernoulli  trials.  If  there  is  no  signal  present, 
since  the  noise  is  determined  by  a  stationary  process,  pf  is  the 
same  for  each  sample  interval  and,  in  terms  of  these  outcomes, 
the  sequence  is  a  series  of  repeated  independent  Bernoulli  trials 
and  therefore  the  expected  number  of  trials  between  false  alarms 
is  1/Pf.  Since  the  time  between  trials  is  St,  the  expected 
number  of  seconds  between  false  alarm  is  St/Pf  or  the  expected 
number  of  false  alarms  per  second  R  is  equal  to  pf/St. 

The  determination  of  x*  depends  on  the  encounter  model 
statistic.  For  both  encounter  models,  when  there  is  a  signal, 

Yl  =  nx  +  alt  y2  -  n2  +  y2,  • • •  >  Ym  *  nm  +  sm  where  s1#  s2, 

• • •  ,  sm  are  signal  values  (voltages)  input  to  a  electrometer’s 
detector.  The  models  imply  that  the  signal  values  s^  =  K-(Es)i 
where  K  is  a  constant  whose  value  is  determined  by  the 
characteristics  of  the  encounter  electrometer  and  where  the 
(Es)i  are  dipole  electric  signal  intensities.  The  models  also 
imply  that  the  noise  values  n^  are  determined  by  a  gaussian 
stochastic  processes  characterized  by  a  standard  deviation  a 
and  that  n^  =  K- (E^) i  +  n|  where  the  (En)i  are  electric  noise 
intensities  and  the  n|  are  electrometer  instrument  noise  values 
(voltages) .  In  the  program,  the  magnitude  of  K  is  1.  Since, 
for  both  models,  p^  depends  only  on  the  ratio  of  signal  energy 
to  noise  energy  for  a  sample  interval,  this  is  a  satisfactory 
choice  for  the  program  if  instrument  noise  is  assumed  to  be 
determined  by  a  process  independent  of  the  electric  noise  process 
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and  to  be  expressed  in  terms  of  an  equivalent  electric  noise 
(E)l  =  ( 1/K)  n[. 

For  both  encounter  models,  the  signal  (the  average  signal 
power)  S  =  (1/m) -E  s£  and  the  noise  (the  expected  value  of  the 
average  noise  power)  N  =  a*  so  that  the  signal-to-noise  ratio 
is  (1/m) -E  s|/a*  where  the  sum  index  i  =  1,  2,  •••  ,  m. 

The  Crosscorrelation  Detector  Statistic:  The  statistic  for  the 
first  encounter  model  is  a  crosscorrelation  detector  statistic 
that  is  defined  by  the  sum 

x  =  E  yi-si 

where  the  summation  index  i  =  1,  2,  • • •  ,  m  and  the  sum  is  over 
the  values  corresponding  to  a  sample  interval.  For  the  first 
encounter  model,  the  characteristics  of  both  the  noise  and  the 
signal  are  required  in  order  to  determine  encounter  detection 
probabilities.  In  particular,  the  signal  values  for  an  encounter 
are  in  the  memory  of  the  detector  prior  to  the  encounter.  For 
the  encounter  conditions  and  a  specified  false  alarm 
probability,  the  statistic  is  optimum  in  the  sense  that  the 
encounter  detection  probability  for  this  statistic  is  at  least 
equal  to  that  for  any  other  statistic.  Because  of  these 
considerations,  encounter  probabilities  based  on  the 
crosscorrelation  statistic  can  be  considered  to  represent  upper 
bounds  on  detection  performance  against  dipole  targets  for 
electrometers  of  the  type  described  by  the  models. 

For  a  sample  interval  without  signal  energy,  x  is  the  value 
of  a  normal  random  variable  with  a  mean  =  0  and  a  variance 


<?X  =  o**S  S£  where  o  is  the  standard  deviation  associated  with 
the  noise  process  and  the  sum  index  1=1,  2,  • • •  ,  m  and  the 
sum  is  over  the  values  corresponding  to  the  sample  interval. 

This  implies  that 

pf  =  1  -  P(x*/ Ox) 

where  P(z)  is  the  standard  normal  cumulative  distribution 
function.  This  relation  is  the  basis  for  determining  the 
threshold  value  x*. 

For  the  sample  interval  with  signal  energy,  x  is  the  value 
of  a  normal  random  variable  with  a  mean  |*x  =  2  si  where  the  sum 
index  i  =  1 ,  2 ,  • • •  ,  m  and  the  sum  is  over  the  values 
corresponding  to  the  sample  interval.  This  implies  that 

Pd  -  1  “  P(v*  "  d*) 

where  v*  =  x*/ax  and  d  =  Z  s^/o*  =  (1/m)  S/N.  This  relation 
is  the  basis  for  determining  encounter  detection  probabilities 
for  the  first  encounter  model.  The  relation  implies  that  for  a 
specified  false  alarm  probability  pf  the  detection  probability 
Pd.  is  an  increasing  function  of  the  signal  to-noise  ratio  S/N. 
The  Square  Law  Detector  Statistic:  The  statistic  for  the  second 
encounter  model  is  a  square  law  (energy)  detector  statistic  that 
is  defined  by  the  sum 

x  =  Z  yj 

where  the  sum  index  i  =  1 ,  2 ,  • • •  ,  m  and  the  sum  is  over  the 
values  corresponding  to  the  sample  interval .  For  the  second 
encounter  model,  only  the  characteristics  of  the  noise  are 
required  to  determine  encounter  detection  probabilities. 
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For  a  sample  interval  without  signal  energy,  x/o*  is  the 


value  of  a  chi-square  random  variable  with  m  degrees  of 
freedom.  This  implies  that 

pf  =  1  -  P (x*/at  |n) 

where  P(x*/a*l®)  is  the  chi-square  cumulative  distribution 
function  for  a  chi-square  random  variable  with  m  degrees  of 
freedom  and  where  o  is  the  standard  deviation  associated  with 
the  noise  process.  This  relation  is  the  basis  for  determining 
the  threshold  value  x*. 

For  the  sample  interval  with  signal  energy,  x/o*  is  the 
value  of  a  noncentral  chi-square  random  variable  with  m 
degrees  of  freedom  and  noncentral  parameter  2  s^/o*  where  the 
sum  index  i  =  1,  2,  • • •  ,  m  and  the  sum  is  over  the  values 
corresponding  to  the  sample  interval.  This  implies  that 

pd  *  1  -  P(xVo’|m,  2  Si/a*) 

where  P(x*/a*|m,  2  s(/a*)  is  the  noncentral  chi-square 
cumulative  distribution  function  for  a  noncentral  chi-square 
random  variable  with  m  degrees  of  freedom  and  noncentral 
parameter  2  sj[/o  =  (l/m)  S/N.  This  relation  is  the  basis  for 
determining  encounter  detection  probabilities  for  the  second 
encounter  model.  The  relation  implies  that  for  a  specified  false 
alarm  probability  pf,  the  detection  probability  pd  is  an 
increasing  function  of  the  signal-to-noise  ratio  S/N.  This  is 
made  more  evident  by  the  following  relation: 

P(x*/a*|m,  2  s£/a* )  =  2{  (a^/j  )  •  exp(-a)  •  P[x*/o*  I  (m  +  2j)]} 
where  the  sum  index  j  =0,  1,  2,  • • •  and  the  parameter 
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a  =  (1/2)  S  s i/o*  =  (m/2) • (S/N)  where  the  sum  index  i  =  1,  2, 
•  •  •  ,  m.  Note  that  the  relation  p<j  >  pf  is  satisfied,  since 
P(x*/at |m)  <  P[x*/0“l(®  +  2  - j ) 3  for  j  =  0,  1,  2,  •  ••  . 
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Appendix  2.  Program  Probability  Calculations 

The  program  evaluates  the  cumulative  and  inverse  cumulative 
distribution  functions  using  approximations  that  are  described  in 
Reference  2.  These  approximations  are  listed  below. 

The  Standard  Normal  Cumulative  Distribution  Function 
Approximation : 

P ( 2 )  =  1  -  st-fbi  +  t-(b2  +  t • (b3  +  t- (b4  +  t • b5) ) ) ) 
where  s  =  (l/(2"jr)  -exp(-z*/2)  and  t  =  1/(1  +  b0z).  And 
where 

b0  =  .2316419,  bL  =  .319381530,  b2  =  -.356563782, 
b4  =  -1.821255978,  b5  =  1.330274429, 
and  z  >  0 .  For  z  <  0,  P(z)  =  1  -  P(|z|). 

The  Inverse  Standard  Normal  Cumulative  Distribution  Function 
Approximation : 

z(P)  =  t  -  (c0  +  t-(cx  +  t  c2)  )/(l  +  t-(d]_  +  t  (d2  +  t  •  d3) )  ) 
where  t  =  (ln(l/Q*))*  and  Q  =  1  -  P.  And  where 
C0  =  2.515517,  cx  =  .802853,  c2  =  .010328, 
dx  =  1.432788,  d2  =  .189269,  d3  =  .001308, 
and  .5  <  P  <  1.  For  0  <  P  <  .5,  z (P)  =  -z (1  -  P) . 

The  Inverse  Chi-Sguare  Cumulative  Distribution  Function 
Approximation : 

v(P|m)  =  m[l  -  2/ (9 ■ m)  +  z • (2/ (9  m) ) i ] 3 
where  P(z)  =  P(v|m) .  In  the  program,  the  inverse  standard 
normal  cumulative  distribution  function  approximation  is  used  to 
determine  z . 
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The  Noncentral  Chi-Square  Cumulative  Distribution  Function 
Approximation : 

P(w|m,E  =  P(z) 

where  z  =  [2-w/(l  +  b)  ] }  -  [2a/(l  +  b)  -  l]i  with 

a  =  m  +  s  s [/a* ,  b  =  (S  Si/o*)/(m  +  E  s^/a*)  and  the  sum  index 
i  *  1,  2,  ••■,m.  In  the  program,  the  standard  normal 
cumulative  distribution  function  approximation  is  used  to 
determine  P(z) . 
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Appendix  3.  The  Electric  Signal 

The  encounter  models  are  defined  by  the  following 
conditions:  A  submarine  electric  anomaly  field  is  an  electric 
dipole  field.  An  electrometer  electric  signal  value  is  the 
magnitude  of  the  component  of  a  dipole  electric  field  in  the 
direction  of  the  electrometer  axis.  The  basis  for  determining 
electric  signal  values  in  the  program  is  an  expression  that  can 
be  developed  as  follows:  In  the  right  handed  rectangular 
coordinate  system  that  is  shown  in  Figure  1,  an  electric  dipole 
is  at  the  origin,  the  xy-plane  is  the  horizontal  plane  at  a 
representative  point  in  an  encounter  region,  the  positive  y-axis 
is  in  the  direction  t  of  the  electrometer's  axis  and  the 
positive  z-axis  is  positive  upward.  In  this  rectangular 
coordinate  system,  the  unit  vector  with  the  orientation  of  the 
electrometer's  axis  can  be  expressed  as  e  =  j-cos  *  +  k-sin  * 
where  #  is  the  electrometer  axis  depression  angle. 

In  a  spherical  coordinate  system  with  the  origin  at  the 
electric  dipole  and  the  polar  axis  in  the  direction  of  the  dipole 
moment,  =  (c-p/r3) • (2 ■ r- cos  ©  +  9- sin  9)  is  the  electric 

field  of  the  dipole  at  a  point  with  spherical  coordinates 
(r,r,9).  In  this  expression  p  is  the  magnitude  of  the  dipole 
moment  and  c  is  a  constant  whose  value  is  determined  by  the 
choice  of  units.  The  electric  signal  for  an  electrometer  that  is 
described  by  the  encounter  models  is  Es  =  e^E^  and  it  can  be 
expressed  in  terms  of  the  rectangular  coordinate  system  as 


» 


« 


« 
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Figure  1.  A  unit  vector  in  the  direction  of  the  electrometer  axis 
and  a  unit  vector  in  the  direction  of  the  dipole  field 
of  a  dipole  at  the  origin  are  shown  at  a  point  that  is 
a  distance  r  from  the  dipole.  The  un^t  vector  h  is 
in  the  direction  of  the  horizontal  component  of  the 
dipole  moment. 


follows:  First,  let  rg  be  a  unit  vector  with  the  orientation 

of  the  spherical  coordinate  system  polar  axis.  Then  the  electric 
dipole  field  =  (c-p/r3) • (3 • r- cos  e  -  r0)  since  the  unit 
vector  r  =  (i-x  +  j-y  +  k-z)/r  where  r  =  (x*  +  y*  +  z’)^  , 

8  =  [(rg  x  r)  x  r]/sin  e  =  [(r*rg)r  -  (r*r) • r0]/sin  8,  and  the 
dot  product  r*r0  =  cos  0.  The  unit  vector  r0  can  be  expressed 
in  the  rectangular  coordinates  by  noting  that  r0  =  p/p  and  then 
expressing  p  in  rectangular  coordinates.  To  do  this,  let  n 
be  the  depression  angle  of  p  from  the  xy-plane  (the  horizontal 
plane)  with  n  positive  downward  and  let  a  be  the  direction  of 
p  relative  to  the  electrometer  axis  direction.  Then  the 
electric  dipole  moment  p  =  p- (h-cos  n  -  k-sin  n)  in  terms  of 
the  unit  vector  k  and  the  unit  vector  h  =  i-sin  a  +  j  •  cos  n 
which  has  the  direction  of  the  horizontal  component  of  p 
relative  to  the  y-axis.  Therefore,  the  unit  vector  in  the  dipole 
direction  r0  =  i  (cos  n-sin  a)  +  j-(cos  fl-cos  a)  -  k-sin  n  and 
Es  =  (c ■ p/r3) • [ (3/r) • cos  0-(ycos  i  -  z-sin  $) 

(cos  i  ■  cos  fi  cos  a  +  sin  i  sin  n)  ] 
where  cos  8  =  (l/r)-(xcos  n-sin  a  +  y-cos  n-cos  a  -  z-sin  n) 
since  cos  0  =  r*r0.  As  can  be  seen  from  this  expression,  for  a 
constant  dipole  moment  magnitude  and  direction  and  a  constant 
electrometer  direction,  the  electrometer  signal  is  a  function  of 
the  rectangular  coordinates  of  the  location  of  the  electrometer 
relative  to  the  dipole.  In  the  encounter  models,  both  of  these 
conditions  are  satisfied. 


Appendix  4.  The  Anderson  Formulation 

In  the  encounter  models,  the  electric  signal  Es  at  a 
sample  point  in  a  straight  line  encounter  can  be  represented  in  a 
form  described  by  Anderson  in  Reference  4.  For  convenience,  the 
Anderson  formulation  is  used  in  the  program  to  determine  values 
for  Es.  It  can  be  developed  as  follows:  In  Figure  2,  the 
primed  rectangular  coordinate  system  is  superimposed  on  the 
rectangular  coordinate  system  of  Figure  1  so  that  the  origin  is 
coincident  with  the  origin  of  that  system.  An  electrometer  is  in 
a  straight  line  encounter  with  an  electric  dipole  that  is  located 
at  the  origin  of  the  combined  system.  The  combined  system  moves 
with  the  electric  dipole  with  the  x'-axis  oriented  so  that  it  is 
parallel  to  and  in  the  direction  of  the  electrometer's  track 
relative  to  the  dipole  and  the  z'-axis  oriented  so  that  it  is 
directed  toward  and  passes  through  the  CPA  on  that  track.  In  the 
primed  coordinate  system,  let  1,  m,  and  n  be  the  direction 
cosines  of  the  dipole  moment  p  and  lj_,  and  n^  be  the 

direction  cosines  of  the  unit  vector  e.  Then,  the  unit  vector 
r0  =  i'l  +  j'n  +  k'  m.  And,  for  points  on  the  relative  track, 
x'  =  s',  y'  =  0  and  z'  =  R  where  R  is  the  slant  range  of  the 
dipole  at  CPA  and  s'  is  the  electrometer's  algebraic  distance 
from  CPA  on  the  relative  track.  (It  is  negative  for  points  before 
CPA  and  positive  for  points  after  CPA.)  This  implies  that  the 
unit  vector  r  =  i' • (s'/r)  +  k' • (R/r) .  From  Appendix  3,  the 
field  is  E(j  =  (c-p/r3)  •  (3  ■  r- cos  0  -  tq)  .  With  r  and  r0 
expressed  in  terms  of  the  primed  unit  vectors  and 
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Figure  2 .  The  dipole  is  at  the  origin  in  the  unprimed  and  primed 
coordinate  systems.  In  the  primed  coordinate  system, 
the  CPA  is  at  (0,  0,  R)  and,  for  a  time  t  during 
the  encounter,  the  electrometer  is  at  [s'(t),  0,  R] . 
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cos  6  =  ttq  =  l-(s'/r-)  +  n-(R/r),  this  becomes 

Ed  =  (c-p/r3) • [ (3/r2) • (Is'  +  n-  R) ■ (i1 -s'  +  k'  R) 

-  (i'-l  +j’-m  +  k* • n) ] . 

Then,  since  e  =  i'  +  j'  m^  +  k'  n^  and  Es  =  e*Ed, 

Es  =  (c  p/r3)  [(2  11!  -  m-mj.  -  n^)-^')* 

+  3-(n-lj_  +  l-nj^-s'-R  +  (2-n-ni  -  l-l^  -  m-mj)  •  R*  ] 
The  quantities 

A2  =  2-1-1^  -  U'jij  -  n-n^ 

Ai  =  3  - (n  - li  +  1 • n^) 

A0  =  2-n  n^  -  l  li  “  Hi -mi, 

are  Anderson  like  coefficients.  With  r=  [(s')*+R*]i  and 
p  =  s'/R,  Es  can  now  be  expressed  as  follows: 

ES(P)  -  (c-p/R3) -2  Aj-Fj(p) 

where  the  Fj(P)  =  pj/(l  +  P*)5/2  for  j  =0,  1,  2  are 
Anderson  like  functions.  This  is  the  Anderson  formulation. 

To  relate  the  Anderson  formulation  for  Es  to  the 
formulation  for  Es  in  Appendix  3,  first  note  that 

1  =  ro*i' 
m  =  r0  •  j  1 

=  r*-.  •  lr  f 


Then  express  r0,  e  and  the  unit  vectors  i',  j'  and  k'  in 
terms  of  the  unit  vectors  i,  j  and  k  and  take  the  indicated 
dot  products.  From  Appendix  3, 

tq  =  i- (cos  fisin  a)  +  j(cos  n-cos  a)  -  k-sin  n 

and 

e  =  j-cos  *  •+•  k  sin  *. 

To  express  i' ,  j'  and  k'  in  terms  of  the  unit  vectors  i,  j 
and  k,  note  that  the  unprimed  coordinate  system  can  be 
transformed  to  the  primed  coordinate  system  by  two  rotations  that 
are  defined  as  follows:  First,  rotate  a  coordinate  system  that 
is  coincident  with  the  unprimed  coordinate  system  about  the 
z-axis  so  that  its  x-axis  is  parallel  to  and  in  the  direction  of 
the  relative  track.  With  positive  values  clockwise  (left  hand 
rule) ,  this  is  a  rotation  through  the  angle  ($  -  n/2)  where  as 
shown  in  Figure  2,  the  direction  of  the  electrometer's  track 
relative  to  the  dipole  is  t  +  <t>.  Next,  rotate  this  system 
about  its  x-axis  through  an  angle  5  with  positive  values 
counterclockwise  (right  hand  rule)  so  that  the  positive  z-axis 
passes  through  the  CPA.  The  angle  S  is  related  to  the 
vertical  separation  z0  of  the  electrometer  and  the  dipole  and 
the  algebraic  encounter  lateral  range  L  that  is  positive  if  the 
dipole  is  to  the  left  of  the  relative  track.  With  these  sign 
definitions:  L  =  Rsin  6  and  z0  =  Rcos  6.  After  the 

rotation,  the  auxiliary  coordinate  system  is  coincident  with  the 
primed  coordinate  system. 
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These  transformations  can  be  described  in  terms  of  matrix 
equations  as  follows:  Let  (xM,y" ,  z")  be  the  coordinates  of  a 
point  in  the  coordinate  system  that  is  coincident  with  the 
auxiliary  coordinate  system  after  the  first  rotation.  Then  the 
transformation  from  the  unprimed  coordinate  system  to  this  double 
primed  coordinate  system  is  described  by  the  matrix  equation 


X" 

sin 

0 

cos 

0 

0 

y" 

= 

-cos 

0 

sin 

0 

0 

z” 

0 

0 

1 

x 

y  • 

z 


And  the  transformation  from  the  double  primed  coordinate  system 
to  the  primed  coordinate  system  is  described  by  the  matrix 
equation 


X' 

10  0 

X" 

y* 

= 

0  cos  6  sin  S 

y" 

z ' 

0  -sin  6  cos  6 

z" 

Taking  the  product  of  the  rotation  matrices  in  the  indicated 
order  yields  the  matrix  equation 

x 

y 

z 


X’ 

sin 

0 

cos 

0 

0 

y’ 

= 

-cos 

6  cos  0 

cos 
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sin  0 

sin 

S 

z  • 

sin 

6  cos  0 

-sin 

5 

sin  0 

cos 

5 

which  defines  the  transformation  from  the  unprimed  to  the  primed 
coordinate  system.  The  unprimed  vector  components  of  the  unit 
vector  i'  can  be  found  by  transforming  the  coordinates  (1,0,0) 


9! 


fi 
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in  the  primed  system  to  their  corresponding  coordinates  in  the 
unprimed  system  with  the  inverse  of  this  matrix  and  then 
repeating  this  process  for  (0,1,0)  and  (0,0,1)  in  order  to  find 
the  unprimed  unit  vector  components  of  j'  and  k' .  However, 
since  the  inverse  transformation  matrix . is  the  transpose . of  this 
matrix,  the  elements  of  the  row  that  corresponds  to  a  primed  unit 
vector  are  the  magnitudes  of  the  unprimed  vectors  that  are  its 
components.  Consequently: 

i»  =  i-sin  <p  +  j  ■  cos  0 

j«  =  -i-cos  S  cos  0  +  j  •  cos  6  sin  0  +  k-sin  6 

k'  =  i-sin  6  cos  0  -  j-sin  6  sin  0  +  k-cos  6. 

Then,  taking  the  dot  products  between  r,  tq  and  these  three 

unit  vectors  as  indicate  above  gives: 

1  =  cos  fi-cos  ( <f>  -  a) 

m  =  cos  6 -cos  fl-sin  (0  -  a)  -  sin  6  sin  n 
n  =  -sin  5-cos  n-sin  (0  -  a)  -  cos  6 -sin  n 

and 

1^  =  cos  *-cos  0 

m^  =  cos  6 • cos  *-sin  0  -  sin  6 ■ sin  # 
n^  =  -sin  6 -cos  *-sin  0  -  cos  S • sin  *. 

These  are  the  relations  that  are  used  in  the  program  to  determine 
values  for  the  Anderson  like  coefficients. 
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Appendix  5.  The  Encounter  Equations  of  Motion 

In  the  double  primed  coordinate  system  that  is  defined  in 
Appendix  4,  the  equations  of  motion  of  an  electrometer  relative 
to  a  submarine  (dipole)  are: 
x*'(t)  =  s' (t) 

y" (t)  =  -l 

2"(t)  =  z0. 

where  L  is  the  algebraic  encounter  lateral  range  that  is 
defined  in  Appendix  4,  z0  is  the  vertical  separation  between  the 
electrometer  and  the  submarine  and  s'(t)  is  the  distance  of  the 
electrometer  from  the  CPA  on  the  relative  track.  With  w  the 
speed  of  the  electrometer  relative  to  the  submarine  and  t  a 
relative  time  parameter,  s'(t)  =  wt  .  These  equations  can  be 
considered  to  be  the  ones  used  in  the  program  to  describe  the 
motion  of  a  electrometer  relative  to  a  submarine.  There,  t  is 
determined  by  t  =  [j  -  (ro-l)/2]-5t  where  the  index  j  =  1,  2, 

• • • ,  m  and  St,  a  time  step,  is  the  time  between  samples.  Note, 
when  t  =  0,  the  electrometer  is  at  the  CPA. 

In  the  coordinate  system  of  Figure  1  in  a  straight  line 
encounter  as  defined  in  the  encounter  models,  the  equations  of 
motion  of  a  electrometer  relative  to  a  submarine  can  be  written 
as  follows: 

x(t)  =  s' (t)  sin  0  +  L-cos  0 
y(t)  =  s' (t) -cos  0  -  L-sin  0 
z(t)  =  z0, 
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since  the  transformation  from  the  double  primed  coordinate  system 
to  the  primed  coordinate  system  is  determined  by  the  matrix 
equation 

x  sin  <t>  -cos  <p  0  x" 

y  =  cos  <fi  sin  <p  0  y" 

z  0  0  1  z " 

x"  =  s' (t)  and  y"  =  -L.  The  above  equations  and  the  expression 

for  Es  in  Appendix  3  could  have  been  used  in  the  program  to 
evaluate  the  electric  signal.  In  particular,  with  these 
equations  of  motion  and  with  the  two  relations  L  =  R-sin  8  and 
z0  =  Rcos  8,  the  expression  for  Hs  in  Appendix  3  can  be 
written  in  terms  of  n,  a,  <J>,  6  and  R  so  that  it  is 
identical  in  appearance  to  the  Anderson  formulation  for  Es  in 
terms  of  these  quantities.  The  definition  of  8  in  Appendix  4 
in  terms  of  a  counterclockwise  rotation  results  in  a  definition 
of  the  algebraic  lateral  range  that  is  consistent  with  some  that 
have  been  used  elsewhere. 

In  the  program,  the  relation  w  =  v  -  u  is  the  basis  for 
determining  the  relative  speed  w.  In  this  relation,  v  is  the 
velocity  of  the  electrometer,  u  is  the  velocity  of  the  submarine 
(dipole)  and  w  is  the  velocity  of  the  electrometer  relative  to 
the  submarine.  This  relation  implies  the  following  equations: 
wx  =  vx  -  ux,  Wy  =  Vy  -  Uy  and  wz  =  vz  -  uz  where  the 
coordinates  x,  y  and  z  refer  to  a  fixed  coordinate  system 
with  the  same  orientation  as  that  of  Figure  1.  In  the  encounters 


0 


of  the  models,  vx  =  v-sin  (a  -  t),  Vy  =  vcos  (o  -  t)  and  vz  = 
where  o  is  the  course  of  the  electrometer  and  v  is  the  speed 
of  the  electrometer.  And,  ux  =  u-sin  (P  -  t) ,  Uy  =  u-cos  (p  -  t) 
and  u2  =  0  where  p  is  the  course  and  u  is  the  speed  of  the 
submarine.  The  angle  <p  and  the  relative  speed  of  the 
electrometer  are  defined  by  wx  =  w-sin  <p  and  wy  *  w-cos  <p. 

In  the  program,  <p  and  w  are  determined  with  a  rectangular  to 
polar  conversion  routine  where  w  =  (  wx  +  Wy)  *  and  <p  is 
determined  by  sin"1(wx/w)  and  cos"1 (Wy/w) . 
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Appendix  6.  The  Submarine  Electric  Dipole 

In  the  encounter  models,  the  magnitude  and  direction  of  a 
submarine's  dipole  moment  p  are  determined  by  first  determining 
its  components  in  the  rectangular  coordinate  system  of  Figure  1. 
In  that  coordinate  system, 

Px  =  PLpsin  (P  -  t)  +  pTp-cos  (p  -  t) 

Py  *  Plpcos  O  ”  t)  -  pTP- sin  (P  -  t) 

Pz  =  “PVP 

where  p  -  t  is  the  submarine's  course  relative  to  the 
electrometer  axis  direction  and  Plp»  Ptp  an(*  PVP  are  t^ie 
longitudinal,  transverse  and  vertical  electric  dipole  moments  of 
the  submarine.  These  relations  are  based  on  the  following  sign 
convention: 

PL  is  positive  when  pL  is  directed  from  stern  to  bow. 

Pt  is  positive  when  pr  is  directed  from  port  to  starboard. 

pv  is  positive  when  py  is  directed  downward. 

The  electric  dipole  moments  are  input  parameters  in  the  program. 
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Appendix  7.  An  Alternative  Encounter  Model 

The  alternative  model  that  is  referred  to  in  Section  III  is 
described  in  more  detail  in  this  appendix.  In  the  model,  the 
detection  range  is  the  slant  range  R  at  the  CPA  for  an 
encounter  with  a  specified  detection  probability  (usually  .5)  is 
defined  by: 

R  =  [cp/Hs]V3 

where  c  and  p  are  defined  in  Appendix  3  and  Es  represents 
a  minimum  detectable  average  electric  signal  that  is  defined  by: 
Es  =  (ORF)  NE 

where  ORF  is  a  signal-to-noise  ratio  called  the  operator 
recognition  factor  and  NE  represents  the  electric  noise. 
Combining  these  two  relations  gives: 

R  =  { c •  p/ [  ( ORF )  NE]}V3. 

The  value  for  ORF  depends  both  on  the  specified  detection 
probability  and  on  a  specified  or  implied  false  alarm 
probability. 

The  Anderson  formulation  is  consistent  with  these  relations 
in  an  approximate  sense  if  the  average  electric  signal  E  is 
defined  as  a  root  mean  square  value  such  that  E  =  (c-p/R3)  k 
where  k  is  an  encounter  parameter  defined  by: 
k  =  {S  [E  Aj-FjOi)]*)* 

with  the  first  sum  index  i  =  1,  2,  • • •  ,  m  and  the  second  sum 
index  j  =  0,  1,  2.  For  a  particular  encounter  geometry,  k  is 
constant  and  this  suggests  that  the  two  encounter  models  could  be 
used  to  determine  an  average  value  for  k  for  an  encounter 
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region  based  on  average  submarine  electric  characteristics. 

Values  for  both  k  and  E  are  generated  by  the  program  and  such 
values  can  give  an  indication  of  the  magnitude  of  the  differences 
in  detection  range  estimates  that  are  based  on  this  model  and 
those  that  are  based  on  either  of  the  other  two  encounter  models. 


34 


Appendix  8.  An  Example  of  the  Program  Output 

The  program  is  designed  to  generate  the  following 
quantities:  encounter  parameter  values,  lateral  range  function 
values  for  the  crosscorrelation  encounter  model  and  for  the 
square  law  encounter  model,  average  electric  signal  values,  slant 
range  at  CPA  values,  encounter  parameter  values  and  electric 
signal  values.  These  values  can  be  saved  as  a  program  data  file 
and/or  they  can  be  printed. 

An  example  of  the  program's  printed  output  is  listed  in 
Table  1,  Table  2  and  Table  3.  Figure  3  is  a  plot  of  the  lateral 
range  function  values  (a  lateral  range  curve)  for  the  square  law 
detector  that  are  listed  in  Table  2.  Figure  4  is  a  plot  of  the 
electric  signal  values  that  are  listed  in  Table  3. 

A  second  example  of  the  program's  printed  output  is  listed 
in  Table  4,  Table  5  and  Table  6.  Figure  5  is  a  plot  of  the 
lateral  range  function  values  for  the  square  law  detector  that 
are  listed  in  Table  5.  Figure  6  is  a  plot  of  the  electric  signal 
values  that  are  listed  in  Table  6. 
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program  file  name 
program  data  file  name 
electric  data  file  name 
processing  data  file  name 
kinematic  data  file  name 


EAD.BAS 

data.ead 

data.ele 

data.prc 

data.kim 


combined  electric,  processing  &  kinematic  data  file  name  data.epk 


electrometer  axis  direction  (decimal  degrees)  335 

electrometer  axis  depression  angle  (decimal  degrees)  0 

encounter  medium  conductivity  (mhos/meter)  4 . 2 

electric  longitudinal  moment  (ampere-meters)  25 

electric  transverse  moment  (ampere-meters)  0 

electric  vertical  moment  (ampere-meters)  0 

sampling  period  (seconds)  10 

integration  time  (seconds)  400 

adjusted  integration  time  (seconds)  410 

number  of  samples  per  encounter  41 

electrometer  course  (decimal  degrees)  0 

electrometer  speed  (knots)  0 

electrometer  altitude  (meters)  0 

target  course  (decimal  degrees)  57 

target  speed  (knots)  10 

target  depth  (meters)  -100 

electrometer  relative  course  (decimal  degrees)  237 

electrometer  relative  speed  (knots)  10 

electrometer-target  vertical  separation  (meters)  -100 

electric  dipole  moment  (ampere-meters)  25 

dipole  moment  azimuth  (decimal  degrees)  57.00003 

dipole  moment  depression  angle  (decimal  degrees)  0 

distance  between  samples  on  the  relative  track  (meters)  51.44445 
false  alarm  rate  (false  alarms  per  hour)  .15 

false  alarm  probability  1.666667E-02 

electric  noise  (volts/meter)  3.5E-08 

maximum  lateral  range  (meters)  500 

lateral  range  step  (meters)  20 

number  of  lateral  range  function  values  51 


Table  1.  An  example  of  an  encounter  parameter  values  printout. 
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data.ead  lateral  range  function  values 


L 

p(cc) 

p<sl) 

E 

R 

k 

meters 

volts/meter 

meters 

-500 

.0359629 

1.49507  IE-02 

5 . 95766E-09 

509.902 

3.221905 

-480 

3 . 869068E-02 

1.506978E-02 

6.713175E-09 

490.306 

3.156311 

-460 

4 . 201478E-02 

1. 522109E-02 

7.58  9571E-09 

470.7441 

3.088742 

-440 

4 . 611728E-02 

.0154154 

8.597962E-09 

451.2206 

3.01912-8 

-420 

5.125207E-02 

1 . 566772E-02 

9 . 758901E-09 

431.7407 

2.947377 

-400 

.0577791 

1 . 599973E-02 

1 . 114912E-08 

412.3106 

2.873376 

-380 

6 . 621807E-02 

1 . 644293E-02 

1. 28209E-08 

392.9377 

2.796981 

-360 

7.733236E-02 

1. 70443  IE-02 

1.492985E-08 

373.6308 

2.718022 

-340 

9 . 226168E-02 

1. 787643E-02 

1.743602E-08 

354.4009 

2.636287 

-320 

.1127298 

1.905402E-02 

2.041433E-08 

335.2611 

2.551514 

-300 

.1413591 

2 . 076598E-02 

2 . 415458E-08 

316.2278 

2.463387 

-280 

.1821076 

2 . 333736E-02 

2.882961E-08 

297.3214 

2.371516 

-260 

.2407389 

2 . 735578E-02 

3 . 454627E-08 

278.5678 

2.275416 

-240 

.3248692 

3 . 394424E— 02 

4.222557E-08 

260 

2.174487 

-220 

.4421915 

4 . 537476E-02 

5 . 17294E-08 

241.6609 

2.067975 

-200 

.5940236 

6 . 647302E-02 

6. 409887E-08 

223.6068 

1.954923 

-180 

.7625715 

.107666 

8 . 065302E-08 

205.9126 

1.834113 

-160 

.9040051 

.1900176 

1. 010635E-07 

188.6796 

1.703969 

-140 

.9785036 

.3451494 

1. 25225E-07 

172.0465 

1.562464 

-120 

.9979309 

.5807582 

1 . 576369E-07 

156.205 

1.406982 

-100 

.9999207 

.8128465 

1 . 953249E-07 

141.4214 

1.234252 

-80 

.9999976 

.9348596 

2 . 313523E-07 

128.0625 

1.04051 

-60 

.9999995 

.9610708 

2 . 582  IE-07 

116.619 

.8225375 

-40 

.9999842 

.8823538 

2 . 39062E-07 

107.7033 

.5807222 

-20 

.9799458 

.3525332 

1 . 489415E-07 

101.9804 

.3276647 

0 

.5137237 

5 . 4 17071E-02 

7.803533E-08 

100 

.1598169 

20 

.9799458 

.3525332 

1. 489415E-07 

101.9804 

.3276647 

40 

.9999842 

.8823538 

2 . 39062E-07 

107.7033 

.5807222 

60 

.9999995 

.9610707 

2 • 5821E-07 

116.619 

.8225374 

80 

.9999976 

.9348597 

2.313523E-07 

128.0625 

1.04051 

100 

.9999207 

.8128465 

1 . 953249E-07 

141.4214 

1.234252 

120 

.9979309 

.5807582 

1. 57637E-07 

156.205 

1.406983 

140 

.9785036 

.3451494 

1.25225E-07 

172.0465 

1.562464 

160 

.9040054 

.1900176 

1.010635E-07 

188.6796 

1.70397 

180 

.7625715 

.107666 

8 . 065302E-08 

205.9126 

1.834113 

200 

.5940236 

6 . 647302E-02 

6. 409888E-08 

223.6068 

1.954923 

220 

.4421915 

4.537467E-02 

5. 172939E-08 

241.6609 

2.067974 

240 

.3248694 

3 . 394424E-02 

4 . 222557E-08 

260 

2.174487 

260 

.2407389 

2 . 735578E-02 

3 . 454627E-08 

278.5678 

2.275416 

280 

.1821077 

2 . 333736E-02 

2 . 882961E-08 

297.3214 

2.371516 

300 

.1413592 

2 . 076598E-02 

2 . 415459E-08 

316.2278 

2.463388 

320 

.1127298 

1 . 905402E-02 

2 . 041433E-08 

335.2611 

2.551514 

340 

9 . 226175E-02 

1 . 787643E-02 

1 . 7  4  3602E-08 

354.4009 

2.636287 

360 

7..733236E-02 

1 . 70443  IE-02 

1 . 492985E-08 

373.6308 

2.718022 

380 

6.62 1807E-02 

1.644293E-02 

1 . 28209E-08 

392.9377 

2.796981 

400 

5 . 777909E-02 

1. 599973E-02 

1 . 114912E-08 

412.3106 

2.873375 

420 

5 . 125207E-02 

1. 566772E-02 

9 . 758901E-09 

431.7407 

2.947377 

440 

4 . 611728E-02 

.0154154 

8 . 597962E-09 

451.2206 

3.019128 

460 

4 . 201478E-02 

1.522109E-02 

7 . 589571E-09 

470.7441 

3.088743 

480 

3.869068E— 02 

1 . 506978E-02 

6.71317  5E-09 

490.306 

3.156311 

500 

.0359629 

1. 495071E-02 

5 . 95766E-09 

509.902 

3 . 221906 

Table  2 


An  example  of  a  lateral  range  function  values  printout 
The  heading  for  the  crosscorrelation  values  is  p(cc) 
and  the  heading  for  the  square  law  values  is  p(sl). 


Figure  3 .  A  plot  of  the  square  law  lateral  range  function  values 
that  are  listed  in  Table  2.  The  horizontal  axis  is 
encounter  horizontal  range  at  CPA  in  meters.  The 
vertical  axis  is  encounter  detection  probability. 
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data.ead  electric  signal  values 
tor  a  lateral  range  of  40  meters 


relative  CPA  distance 
s  in  meters 

-1028.889 

-977.4446 

-926.0001 

-874.5557 

-823.1112 

-771.6667 

-720.2223 

-668.7778 

-617.3334 

—565 . 889 

-514.4445 

-463.0001 

-411.5556 

-360.1111 

-308.6667 

-257.2223 

-205.7778 

-154.3333 

-102.8889 

-51.44445 

0 

51.44445 

102.8889 

154.3333 

205.7778 

257.2223 

308.6667 

360.1111 

411.5556 

463.0001 

514.4445 

565.889 

617.3334 

668.7778 

720.2223 

771.6667 

823.1112 

874.5557 

926.0001 

977.4446 

1028.839 


electric  signal 
Es  in  volts/aeter 

1 . 660273E-10 
1 . 959909E-10 
2 . 335016E-10 
2 . 8 10548E-10 
3 . 42 1893E-10 
4.220248E-10 
5 . 28136E-10 
6 . 720149E-10 
8 . 715932E-10 
1.155758E-09 
1 . 572761E-09 
2 . 206581E-09 
3 . 210406E-09 
4 . 878809E-09 
7.B1198E-09 
1 . 330549E-08 
2.426627E-08 
4 . 685768E— 06 
8.75983E-08 
1 . 023685E-07 
-5 . 276525E-08 
-1. 366934E-07 
-7.038444E-08 
-2 . 673259E-08 
-1. 000855E-08 
-3 . 866346E-09 
—1. 493584E— 09 
-5 . 247395E-10 
-1. 147233E-10 
5 . 922923E-11 
1 . 2893  36E-10 
1 . 514067E-10 
1 . 524577E-10 
1.44157E-10 
1.32254E-10 
1 . 194879E-10 
1 . 071295E-10 
9 . 572235E-11 
8 . 545079E-11 
7 . 632433E-11 
6. 827254 E- 11 


Table  3.  An  example  of  an  electric  signal  values  printout. 
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-1028  869 


-514  4445 


0 


514.4445 


1028  889 


Figure  4 


.  A  plot  of  the  electric  signal  values  listed  in  Table  3 
The  horizontal  axis  is  relative  CPA  distance  in  meters 
The  vertical  axis  is  electric  signal  in  nanovolts  per 


meter. 


program  file  name  EAD.BAS 

program  data  file  name  datal.ead 

electric  data  file  name  data.ele 

processing  data  file  name  data.prc 

kinematic  data  file  name  datal.kim 

combined  electric,  processing  t  kinematic  data  file  name  datai.epk 

electrometer  axis  direction  (decimal  degrees)  335 

electrometer  axis  depression  angle  (decimal  degrees)  0 

encounter  medium  conductivity  (mhos/meter)  4.2 

electric  longitudinal  moment  (ampere-meters)  25 

electric  transverse  moment  (ampere-meters)  0 

electric  vertical  moment  (ampere-meters)  0 

sampling  period  (seconds)  10 

integration  time  (seconds)  400 

adjusted  integration  time  (seconds)  410 

number  of  samples  per  encounter  41 

electrometer  course  (decimal  degrees)  0 

electrometer  speed  (knots)  0 

electrometer  altitude  (meters)  0 

target  course  (decimal  degrees)  170 

target  speed  (knots)  10 

target  depth  (meters)  -100 

electrometer  relative  course  (decimal  degrees)  350 

electrometer  relative  speed  (knots)  10 

electrometer-target  vertical  separation  (meters)  -100 

electric  dipole  moment  (ampere-meters)  25 

dipole  moment  azimuth  (decimal  degrees)  170 

dipole  moment  depression  angle  (decimal  degrees)  0 

distance  between  samples  on  the  relative  track  (meters)  51.44444 
false  alarm  rate  (false  alarms  per  hour)  .15 

false  alarm  probability  1.666667E-02 

electric  noise  (volts/meter)  3.5E-08 

maximum  lateral  range  (meters)  500 

lateral  range  step  (meters)  20 

number  of  lateral  range  function  values  51 


m 


9 


Table  4. 


fl 

An  example  of  an  encounter  parameter  values  printout. 
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datal.ead  lateral  range  function  values 


L 

p(cc) 

P(sl) 

E 

R 

k 

meters 

volts/meter 

meters 

-500 

3 . 110245E-02 

1 . 475347E-02 

4 . 65646E-09 

509.9019 

2.582769 

-480 

3 . 308339E-02 

1.483123E-02 

5.235783E-09 

490.306 

2.536152 

-460 

3 . 548338E-02 

1 . 493034E-02 

5 . 903325E-09 

470.7441 

2.488133 

-440 

3 . 842717E-02 

1.505803E-02 

6. 703113E-09 

451.2206 

2.438668 

-420 

4.208816E-02 

1 . S22449E-02 

7 . 637841E-09 

431.7407 

2.387708 

-400 

4 . 671227E— 02 

1 . 544411E-02 

8 . 748915E-09 

412.3106 

2.335211 

-380 

5.265528E-02 

1. 573816E-02 

1. 009386E-08 

392.9377 

2.28113 

-360 

6. 044344E-02 

1 . 613816E— 02 

1. 168881E-08 

373.6308 

2.225419 

-340 

.0708726 

1. 669256E-02 

1 . 3 6 640 IE— 08 

354.4009 

2.168032 

-320 

8.517325E-02 

1 . 747827E-02 

1. 606643E-08 

335.2611 

2.108919 

-300 

.1052861 

1 . 862 119E-02 

1.899775E-08 

316.2278 

2.048029 

-280 

.1343149 

2 . 033773E-02 

2 . 283886E-08 

297.3214 

1.985313 

-260 

.1772254 

2 . 302039E-02 

2 . 769995E-08 

278.5678 

1.92072 

-240 

.2417317 

2.742749E-02 

3 . 402688E-08 

260 

1.854206 

-220 

.3386843 

.0351379 

4 . 239666E-08 

241.6609 

1.785746 

-200 

.4792239 

4 . 971678E-02 

5 . 322737E-08 

223.6068 

1.715341 

-180 

.662185 

7 . 983379E-02 

6 . 817646E-08 

205.9126 

1.643053 

-160 

.8491701 

.1472224 

8 . 849707E-08 

188.6796 

1.569052 

-140 

.9671919 

.300248 

1. 156304E-07 

172.0465 

1.4937 

-120 

.9981669 

.5916718 

1 . 542097E-07 

156.205 

1.417713 

-100 

.9999912 

.9016253 

2.073466E-07 

141.4214 

1.34246 

-80 

1 

.9966924 

2.75907  4E-07 

128.0625 

1.270434 

-60 

1 

.9999966 

3 . 585535E-07 

116.619 

1.205756 

-40 

1 

1 

4 .46597  6E-07 

107.7033 

1.154001 

-20 

1 

1 

5. 20507E-07 

101.9804 

1.120663 

0 

1 

1 

5 . 416029E-07 

100 

1.109207 

20 

1 

1 

5.20507E-07 

101.9804 

1.120663 

40 

1 

1 

4.465976E-07 

107.7033 

1.154001 

60 

1 

.9999966 

3.585535E-07 

116.619 

1.205757 

80 

1 

.9966924 

2.759074E-07 

128.0625 

1.270434 

100 

.9999912 

.9016253 

2 . 073466E-07 

141.4214 

1.34246 

120 

.9981669 

.5916718 

1 . 542097E-07 

156.205 

1.417713 

140 

.9671919 

.300248 

1. 156304E-07 

172.0465 

1.493699 

160 

.8491701 

.1472224 

8.849707E-08 

188.6796 

1.569052 

180 

.662185 

7 . 983379E-02 

6.817645E-08 

205.9126 

1.643053 

200 

.4792239 

4.971678E-02 

5.322737E-08 

223.6068 

1.715341 

220 

.3386843 

.0351379 

4 . 239666E-08 

241.6609 

1.785746 

240 

.2417317 

2 . 742749E-02 

3 . 402688E-08 

260 

1.854206 

260 

.1772254 

2 . 302039E-02 

2 . 769995E-08 

278.5678 

1.920719 

280 

.1343149 

2 . 033773E-02 

2 . 283886E-08 

297.3214 

1.985313 

300 

.1052861 

1 . 862119E-02 

1. 899775E-08 

316.2278 

2.04803 

320 

8.517325E-02 

1 . 747827E-02 

1. 606643E— 08 

335.2611 

2.108919 

340 

7 . 087261E-02 

1 . 669256E-02 

1.366401E-08 

354. 400S 

2.168032 

360 

6 . 044344E-02 

1. 613816E-02 

1 . 16888  IE-08 

373.6308 

2.225419 

380 

5 . 265528E-02 

1 . 573816E-02 

1 . 009386E-08 

392.9377 

2.28113 

400 

4 . 671226E-02 

1 . 5444 11E-02 

8.748914E-09 

412.3106 

2.335211 

420 

4.208816E-02 

1 . 522449E-02 

7 . 637841E-09 

431.7407 

2.387708 

440 

3.842716E-02 

1 . 505803E-02 

6 . 703 113E-09 

451.2206 

2.438667 

460 

3 . 548338E-02 

1 . 493034E-02 

5.903325E-09 

470.7441 

2.488133 

480 

3 . 308338E-02 

1 . 483 123E-02 

5 . 235783E-09 

490.306 

2.536152 

500 

3 . 110243E— 02 

1 . 475347E-02 

4 . 65646E-09 

509.9019 

2.582769 

Table  5.  An  example  of  a  lateral  range  function  values 

printout.  The  heading  for  the  crosscorrelation  values 
is  p(cc)  and  the  heading  for  the  square  law  values  is 
P(sl) . 
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Lateral  Ronge 


I 


A  plot  of  the  square  law  lateral  range  function  values 
that  are  listed  in  Table  5.  The  horizontal  axis  is 
encounter  horizontal  range  at  CPA  in  meters.  The 
vertical  axis  is  encounter  detection  probability. 


datal . ead  electric  signal  values 
for  a  lateral  range  of  40  meters 


relative  CPA  distance 
s  in  meters 

-1028.889 

-977.4443 

-925.9999 

-874.5555 

-823.111 

-771.6666 

-720.2222 

-668.7777 

-617.3333 

-565.8889 

-514.4444 

-462.9999 

-411.5555 

-360.1111 

-308.6666 

-257.2222 

-205.7778 

-154.3333 

-102.8889 

-51.44444 

0 

51.44444 

102.8889 

154.3333 
205.7778 

257.2222 

308.6666 
360.1111 

411.5555 

462.9999 
514.4444 

565.8889 

617.3333 
668.7777 

720.2222 

771.6666 
823.111 

874.5555 

925.9999 
977.4443 
1028.889 


electric  signal 
Es  in  volts/meter 

-8 . 258498E-10 
-9 . 606357E-10 
-1.126187E-09 
— 1. 331725E-09 
-1. 589959E-09 
-1. 918711E-09 
-2 . 34344 2E-09 
-2 . 901288E-09 
-3.647673E-09 
-4 . 667418E-09 
-6 . 093956E-09 
-8 . 143509E-09 
-1. 117728E-08 
-1.581578E-08 
-2 . 314237E-08 
-3 . 500006E-08 
-5. 39567E-08 
-7 . 945557E-08 
—8 . 038128E-08 
8 . 787511E-08 
3 . 662163E-07 
1.503571E-07 
-3 . 909039E-08 
—6 . 022178E-08 
-4 . 499852E-08 
-3 . 051197E-08 
—2 . 071024E— 08 
— 1 . 44 03 5E— 08 
-1 . 03082 IE-08 
-7 . 582271E-09 
-5 . 716593E-09 
-4 . 404917E-09 
-3 . 459717E-09 
-2.763325E-09 
-2 . 239973E-09 
-1.839639E-09 
-1 . 52852  3E-09 
-1 . 283  285E-09 
-1 . 087492E-09 
-9. 293593E-10 
-8 . 003 003 E- 10 


Table  6.  An  example  of  an  electric  signal  values  printout. 
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.889 


—51  4.44-45 


0 


51 4.4445 


1028  889 


A  plot  of  electric  signal  values  listed  in  Table  6. 
The  horizontal  axis  is  relative  CPA  distance  in 
meters.  The  vertical  axis  is  electric  signal  in 
nanovolts  per  meter. 


Appendix  9.  A  Program  Listing 

10  CLS  :  CLEAR 

20  FRIOT  "Electric  Anomaly  Detection  (EAD)  Lateral  Range  Function  Program" 

30  DIM  X0(70) ,  R0(70),  E(70) ,  ES(70,  150),  K(70) ,  PDCC(70) ,  PDSL(70) 

40  PI  =  3.141592654#:  CON  =  2  *  PI  /  360:  KDU  =  1852  /  3600:  N$  =  "EAD. BAS" 

50  Q0  =  .2316419:  Q1  =  .31938153#:  Q2  =  -.356563782#:  Q3  =  1.781477937# 

60  Q4  =  -1.821255978#:  Q5  =  1.330274429# 

70  II  =  2.515517:  12  =  .802853:  13  =  .010328:  14  =  1.432788:  15  =  .189269:  16 
=  .001308 
80  PRINT 

90  INPUT  "generate  data  or  print  a  program  data  file  (g/p) " ;  A$ 

100  IF  A$  =  "G"  OR  A$  =  "g"  THEN  GOTO  120 

110  IF  A$  =  "P"  OR  A$  =  "p"  THEN  GOTO  2700  ELSE  GOTO  90 

120  PRINT  :  A$  =  "a" 

130  INPUT  "electric,  processing  &  kinematic  data  entry  by  combined  file 
(y/n)";  A$ 

140  IF  A$  -  "Y"  OR  A$  =  "y"  THEN  GOTO  1740 

150  IF  A$  =  "N"  OR  A$  =  "n"  THEN  GOTO  160  ELSE  GOTO  130 

160  INPUT  "electric  data  entry  by  file/keyboard  (f/k)";  A$ 

170  IF  A$  =  "F"  OR  A$  =  "f"  THEN  GOTO  430 

180  IF  A$  =  "K"  CR  A$  =  "k"  THEN  GOTO  190  ELSE  GOTO  160 

190  INPUT  "electrometer  axis  direction  in  decimal  degrees";  EAD 

200  INPUT  "electrometer  axis  depression  in  decimal  degrees  (down  +)";  DIP 

210  A$  =  "a" 

220  INPUT  "input  encounter  medium  conductivity  (y/n) " ;  A$ 

230  IF  A$  =  "Y"  OR  A$  =  "y"  THEN  GOTO  250 

240  IF  A$  «  "N"  OR  A$  =  "n"  THEN  GOTO  270  ELSE  GOTO  220 

250  INPUT  "encounter  medium  conductivity  in  mobs";  OOND 

260  GOTO  280 

270  OOND  =  4 

280  ELM  =  0:  ELM  =  0:  EVM  =  0 

290  INPUT  "longitudinal  electric  moment  in  anpere-meters  (stem-to-bcw  +)"; 
ELM 

300  INPUT  "transverse  electric  moment  in  anpere-meters  ( port-to-starboard  +)"; 
ETM 

310  INPUT  "vertical  electric  moment  in  anpere-meters  (downward  +)";  EVM 
320  A$  =  "a" 

330  INPUT  "generate  an  electric  data  file  (y/n)";  A$ 

340  IF  A$  =  "Y"  OR  A$  =  "y"  THEN  GOTO  360 

350  IF  A$  =  "N"  OR  A$  =  "n"  THEN  GOTO  490  ELSE  GOTO  330 

360  INPUT  "electric  data  file  name";  E$ 

370  ON  ERROR  GOTO  380:  GOTO  390 

380  RESUME  360 

390  OPEN  "O",  #1,  E$ 

400  WRITE  #1,  EAD,  DIP,  OOND,  ELM,  ELM,  EVM 
410  CLOSE 
420  GOTO  490 

430  INPUT  "electric  data  file  name";  E$ 

440  ON  ERROR  GOTO  450:  GOTO  460 

450  RESUME  430 

460  OPEN  "I",  #1,  E$ 

470  INPUT  #1,  EAD,  DIP,  OOND,  ELM,  ETM,  EVM 
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480  CLOSE 
490  A$  =  "a" 

500  INPUT  "processing  data  entry  by  file/keyboard  (f/k )";  A$ 

510  IF  A$  =  "F"  OR  A$  =  "f"  THEN  GOTO  1090 

520  IF  A$  =  "K"  OR  A$  =  "k"  THEN  GOTO  530  ELSE  GOTO  500 

530  A$  *  "a" 

540  INPUT  "input  sampling  period  (y/n)";  A$ 

550  IF  A$  =  "Y"  OR  A$  =  "y"  THEN  GOTO  570 

560  IF  A$  =  "N"  OR  A$  =  "n"  THEN  GOTO  600  ELSE  GOTO  540 

570  INPUT  "sanpling  period  in  seconds";  DT 

580  IF  DT  <=  0  THEN  FRLNT  ;  PRINT  "nast  be  greater  than  zero":  PRINT  :  GOTO 
570 

590  GOTO  730 
600  A$  =  "a" 

610  INPUT  "input  maximum  electric  signal  frequency  (y/n)"?  A$ 

620  IF  A$  =  "Y"  OR  A$  =  "y"  THEN  GOTO  640 

630  IF  AS  «  "N"  OR  A$  =  "n"  THEN  GOTO  670  ELSE  GOTO  610 

640  INPUT  "maximum  electric  signed  frequency  in  Hertz"?  MAXF 

650  IF  MAXF  <=  o  THEN  PRINT  :  PRINT  "must  be  greater  than  zero":  PRINT  :  GOTO 

640 

660  GOTO  720 

670  INPUT  "minimum  target  slant  range  at  CPA  in  meters";  MINRO 

680  IF  MINRO  <=  0  THEN  PRINT  :  PRINT  "must  be  greater  than  zero":  PRINT  :  GOTO 

670 

690  INPUT  "maximum  electrometer  relative  speed  in  knots";  MAXVEK 
700  MAXVE  =  MAXVEK  *  KDN:  MAXF  =  2  *  MAXVE  /  MINRO 

710  IF  MAXVM  <=  0  THEN  PRINT  :  PRINT  "must  be  greater  than  zero":  PRINT  :  GOTO 
690 

720  DT  -  1  /  (2  *  MAXF) :  REM  low  pass  filter,  Nyquist  sanpling  rate 
730  A$  =  "a" 

740  INPUT  "input  integration  time  (y/n)";  A$ 

750  IF  A$  =  "Y"  OR  A$  =  "y"  THEN  GOTO  770 

760  IF  A$  =  "N"  OR  A$  =  "n"  THEN  GOTO  850  ELSE  GOTO  740 

770  INPUT  "integration  time  in  seconds";  IT 

780  IF  rr  >=  DT  THEN  GOTO  810 

790  PRINT  :  PRINT  "IT  =  "  +  STR$(IT)  +  "  seconds  -  minimum  =  "  +  STR$  (DT)  +  " 
seconds" :  PRINT 
800  GOTO  770 

810  NS  =  2  *  INT(IT  /  DT  /  2)  +1:  REM  adj  number  of  sanples  per  integration 
time 

820  IF  NS  <=  151  THEN  GOTO  980 

830  PRINT  :  PRINT  "IT  =  "  +  STR$(IT)  +  "  seconds  -  maximum  =  "  +  STR$(150  * 
DT)  +  "  seconds":  PRINT 
840  GOTO  770 

850  INPUT  "maximum  target  slant  range  at  CPA  in  meters";  MAXRO 

860  IF  MAXRO  <=  0  THEN  PRINT  :  PRINT  "must  be  greater  than  zero":  PRINT  :  GOTO 

850 

870  INPUT  "minimum  electrometer  relative  speed  in  knots";  MINVEK 
880  MINVE  =  MINVEK  *  KDN 

890  IF  MINVE  <=  0  THEN  PRINT  :  PRINT  "must  be  greater  than  zero":  PRINT  :  GOTO 
870 

900  IT  =  2  *  MAXRO  /  MINVE 
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910  IF  IT  >=  DT  THEN  GOTO  940 

920  PRINT  :  FEINT  "it  =  "  +  STR$(IT)  +  "  seconds  -  minimum  =  "  +  STR$(DT)  +  " 
seconds":  PRINT 
930  GOTO  740 

940  NS  =  2  *  INT(IT  /  DT  /  2)  +1:  REM  adjusted  number  of  samples  per 

integration  time 

950  IF  NS  <=  151  THEN  GOTO  980 

960  HUNT  :  FEINT  "IT  =  "  +  STE$(IT)  +  "  seconds  -  minimum  =  "  +  STR$(150  * 
DT)  +  "  seconds":  FEINT 
970  GOTO  740 
980  A$  =  "a" 

990  INFUT  "generate  a  processing  data  file  (y/n)";  AS 
1000  IF  A$  =  "Y"  OR  A$  =  "y"  THEN  GOTO  1020 

1010  IF  A$  =  ,rN"  OR  A$  =  "n"  THEN  GOTO  1150  ELSE  GOTO  990 

1020  INPUT  "processing  data  file  name";  P$ 

1030  ON  ERROR  GOTO  1040:  GOTO  1050 
1040  RESUME  1020 
1050  OPEN  "O",  #1,  P$ 

1060  WRITE  #1,  DT,  IT,  NS 

1070  CIOSE 
1080  GOTO  1150 

1090  INFUT  "processing  data  file  name";  P$ 

1100  ON  ERROR  GOTO  1110:  GOTO  1120 
1110  RESUME  1090 
1120  OPEN  "I",  #1,  P$ 

1130  INPUT  #1,  DT,  IT,  NS 

1140  dOSE 
1150  A$  =  "a" 

1160  INPUT  "kinematic  data  entry  by  file/keyboard  (fA)";  A$ 

1170  IF  A$  a  "F"  OR  A$  =  "f"  THEN  GOTO  1440 

1180  IF  A$  =  "K"  OR  AS  =  "k"  THEN  GOTO  1190  EISE  GOTO  1150 

1190  INFUT  "electrometer  course  in  decimal  degrees  (0  if  stationary)";  CE 

1200  CER  -  CE  *  OON:  REM  electrometer  course  in  radians 

1210  INFUT  "eiecrcmeter  speed  in  knots";  VEK 

1220  INFUT  "electrometer  altitude  in  meters  (below  0  is  -) " ;  AE 

1230  INFUT  "target  course  in  decimal  degrees";  CT 

1240  CTR  =  CT  *  OCN:  REM  target  course  in  radians 

1250  INFUT  "target  speed  in  knots";  VTK 

1260  INFUT  "target  depth  in  meters  (above  0  is  -)";  AT 

1270  WXK  =  VEK  *  SIN  (CER)  -  VTK  *  SIN  (CTR) 

1280  WYK  =  VEK  *  006  (CER)  -  VTK  *  COS  (CTR) 

1290  Z  =  AE  +  AT:  REM  vertical  separation  (-  for  electrometer  below  target) 
1300  X  =  WXK:  Y  =  WYK:  GOSUB  3680 

1310  CR  =  T:  C  =  CR  /  OCN:  REM  electrometer  relative  course 
1320  WK  =  R:  REM  electrometer  relative  speed 
1330  A$  =  "a" 

1340  INFUT  "generate  a  kinematic  data  file  (y/n)";  AS 
1350  IF  A$  =  "Y"  CR  A$  =  "y"  THEN  GOTO  1370 
1360  IF  A$  =  "N"  OR  A$  =  "n"  THEN  GOTO  1510  EISE  GOTO  1340 
1370  INFUT  "kinematic  data  file  name";  K$ 

1380  ON  ERROR  GOTO  1390:  GOTO  1400 
1390  RESUME  1370 
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1400  OPEN  "0",  #1,  K$ 

1410  WRITE  #1,  CE,  VEK,  AE,  CT,  VTK,  AT,  Z,  C,  WK 
1420  dOSE 
1430  GOTO  1510 

1440  INPUT  "kinematic  data  file  name";  K$ 

1450  ON  ERROR  GOTO  1460:  GOTO  1470 
1460  RESEME  1440 
1470  OFIN  "I",  #1,  K$ 

1480  INPUT  #1,  CE,  VEK,  AE,  CT,  VTK,  AT,  Z,  C,  WK 
1490  CLOSE 

1500  CIR  =  CT  *  CON:  CR  =  C  *  CON 
1510  EAER  =  EAD  *  OCN 

1520  EMX  =  ELM  *  SIN(Cm  -  EAER)  +  ETM  *  OOS(CER  -  EAER) 

1530  EMY  =  ELM  *  COS  (CIR  -  EAER)  -  ETM  *  SIN  (CIR  -  EAER) 

1540  EMZ  =  -EVM 
1550  X  =  EMX:  Y  =  EMY 
1560  GQSUB  3680 

1570  CMLR  =  T:  REM  dipole  azimuth  relative  to  the  electrometer  axis  direction 
1580  X  =  -EMZ:  Y  =  R 
1590  GQSUB  3680 

1600  EM  =  R:  CMR  =  T:  REM  dipole  depression  angle 
1610  CML  =  CMLR  /  CON:  CM  *  CMR  /  CCN 
1620  A$  =  "a" 

1630  INEUT  "generate  a  combined  electric,  processing  &  kinematic  data  file 
(y/n)";  A$ 

1640  IF  A$  =  "Y"  OR  A$  =  "y"  THEN  GOTO  1660 

1650  IF  A$  =  "N"  OR  A$  =  "n"  THEN  GOTO  1820  ELSE  GOTO  1630 

1660  INPUT  "combined  electric,  processing  &  kinematic  data  file  name";  C$ 

1670  (XL  ERROR  GOTO  1680:  GOTO  1690 
1680  RESUME  1660 
1690  OPEN  "O",  #1,  C$ 

1700  WRITE  #1,  EAD,  DIP,  OOND,  ELM,  ELM,  EVM,  DT,  TT 

1710  WRITE  #1,  NS,  CE,  VEK,  AE,  CT,  VTK,  AT,  Z,  C,  WK,  EM,  CML,  CM,  E$,  P$,  K$ 
1720  CLOSE 
1730  GOTO  1820 

1740  INPUT  "combined  electric,  processing  &  kinematic  data  file  name";  C$ 

1750  ON  ERROR  GOTO  1760:  GOTO  1770 
1760  RESUME  1740 
1770  OPEN  "I",  #1,  C$ 

1780  INPUT  #1,  EAD,  DIP,  COND,  ELM,  ETM,  EVM,  DT,  IT 

1790  INPUT  #1,  NS,  CE,  VEK,  AE,  CT,  VTK,  AT,  Z,  C,  WK,  EM,  CML,  CM,  E$,  P$,  K$ 
1800  CLOSE 

1810  EAER  =  EAD  *  OCN:  CR  =  C  *  CON:  CMLR  =  CML  *  OCN:  CMR  =  CM  *  OCN 
1820  INEUT  "required  false  alarm  rate  in  false  alarms  per  hour";  FAR 
1830  PF  =  EAR  *  IT  /  3600:  REM  false  alarm  probability 
1840  Y  =  PF:  IF  PF  >  .5  THEN  Y  =  1  -  Y:  REM  inverse  normal  approximation 
1850  Y  =  SQR(LOG(l  /  Y  /  Y) ) 

1860  Y  =  Y  -  (II  +  Y  *  (12  +  13  *  Y) )  /  (1  +  Y  *  (14  +  Y  *  (15  +  16  *  Y) ) ) 

1870  IF  PF  <  .5  THEN  Y  =  - Y 
1880  ZP  =  -Y 

1890  CHI  =  NS  *  (1-2/9/NS  +  ZP*  SQR(2  /  9  /  NS) )  *  3:  REM  inverse 
chi-square  approximation 


] 
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1900  INPUT  "electric  noise  in  volts/meter" ;  SIG 
1910  INPUT  "maximum  lateral  range  in  meters";  LRM 
1920  INPUT  "lateral  range  step  in  meters";  ST 
1930  IF  ST  <=  IFM  THEN  GOTO  1950 

1940  PRINT  :  FRENT  "maximum  step  is  "  +  S7R$(LfM)  +  "  meters":  PRINT  :  GOTO 
1920 

1950  NC  =  2  *  3NT(IFM  /  ST)  +1:  REM  number  of  lateral  range  function  values 
1960  IF  NC  <=  71  THEN  GOTO  1980 

1970  PRINT  :  PRINT  "minimum  step  is  "  +  STR$(IPM  /  35)  +  "  meters":  PRINT  : 
GCTO  1920 

1980  AIT  =  DT  *  NS:  REM  adjusted  integration  time  in  seconds 

1990  W  =  WK  *  HON:  REM  electrometer  relative  speed  in  meters/ second 

2000  DS  =  W  *  DT:  REM  distance  between  samples  on  the  relative  track  in  meters 

2010  XO  «  -(NC  -  1)  /  2  *  ST 

2020  FOR  I  =  0  TO  NC  -  1 

2030  X0(I)  =  XO 

2040  X  =  XO:  Y  =  Z 

2050  GQSUB  3680 

2060  RO  =  R:  R0(I)  =  R:  REM  target  slant  range  at  CPA  in  meters 
2070  DELR  =  T:  REM  target  depression  angle  canplement  at  CPA  in  radians 
2080  IF  RO  =  0  THEN  GOTO  2430:  REM  zero  lateral  range  and  vertical  separation 
2090  EMC  =  4  *  PT  *  OOND 

2100  EMF  =  EM  /  EMC  /  RO  *  3:  REM  dipole  mcrrent  factor 

2110  DIFR  =  DIP  *  OON:  REM  dipole  depression  angle  in  radians 

2120  KR  =  CR  -  EADR:  REM  course  referenced  to  the  electrometer  axis  direction 

2130  L  =  GQS(CMR)  *  OOS(KR  -  CMLR) 

2140  M  =  CDS  (DEER)  *  CDS  (OMR)  *  SIN(KR  -  CMIR)  -  SIN  (DELR)  *  SIN(CMR) 

2150  N  =  -SIN (DELR)  *  COS(CMR)  *  SIN(KR  -  CMLR)  -  COS  (DELR)  *  SIN(CMR) 

2160  IA  =  COS(DIPR)  *  COS(KR) 

2170  Ml  *  COS  (DEAR)  *  OOS(DIFR)  *  SIN(KR)  -  SIN  (DELR)  *  SIN  (DIFR) 

2180  N1  =  -SIN (DELR)  *  OOS(DIPR)  *  SIN(KR)  -  CDS (DELR)  *  SIN  (DIFR) 

2190  A2  =  2*L*L1-M*M1-N*  Nl:  REM  anderson  like  function  coefficient 
2200  A1  =  3  *  (N  *  LI  +  L  *  Nl) :  REM  anderson  like  function  coefficient 
2210  AO  =  2  *N*N1-L*L1-M*  Ml:  REM  anderson  like  function  coefficient 
2220  SUM  =  0:  EMAX  =  0:  EMIN  =  0 
2230  FOR  J  =  0  TO  NS  -  1 

2240  S  =  (J  -  (NS  -  1)  /  2)  *  DS:  BA  =  S  /  RO:  REM  anderson  like  function 
argument 

2250  AF  —  1  /  (1  +  BA*  BA)  *  2.5:  REM  anderson  like  function  factor 
2260  ESF  =  (A2  *  BA  *  BA  +  A1  *  BA  +  AO)  *  AF:  REM  electric  signal  factor 
2270  ES(I,  J)  =  EMF  *  ESF:  REM  electric  signal  value 
2280  IF  ES(I,  J)  >  EMAX  THEN  EMAX  =  ES(I,  J) 

2290  IF  ES(I,  J)  <  EMIN  THEN  EMIN  =  ES(I,  J) 

2300  SUM  =  SUM  +  ESF  *  ESF 
2310  NEXT  J 

2320  E(I)  =  EMAX  -  EMIN 
2330  K(I)  =  SQR(SUM) 

2340  W  =  -ZP  +  EMF  *  SQR(SUM)  /  SIG 

2350  LAM  =  EMF  *  EMF  *  SUM  /  (SIG  *  SIG) :  AA  =  NS  +  LAM:  BB  =  1  +  LAM  /  (NS  + 
LAM) 

2360  ZN  =  -SQR(2  *  CHI  /  BB)  +  SQR(2  *  AA  /  BB  -  1) :  XI  =  W 
2370  GOSUB  3740 
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2380  IF  Yl  >  1  THEN  Y1  =  1 

2390  PDCC(I)  =  Yl:  XI  -  ZN 

2400  GQSUB  3740 

2410  IF  Yl  >  1  THEN  Yl  =  1 

2420  PDSL(I)  =  Yl 

2430  X0  =  X0  +  ST 

2440  NEXT  I 

2450  C  =  C  /  360:  C  =  (C  -  INT(C))  *  360 
2460  IF  C  <  0  THEN  C  =  360  +  C 

2470  CMLN  =  (OML  +  EAD)  /  360:  CMLN  =  (CMLN  -  IOT(CMIN))  *  360 
2480  IF  CMLN  <  0  THEN  CMLN  =  360  +  CMLN 
2490  PRINT  :  A$  =  "a" 

2500  INPUT  "generate  a  program  data  file  (y/n)";  A$ 

2510  IF  A$  =  "Y"  OR  A$  =  "y"  THEN  GOTO  2530 

2520  IF  A$  =  "N"  CR  AS  =  "n"  THEN  GOTO  2860  ELSE  GOTO  2500 

2530  INPUT  "program  data  file  name";  D$ 

2540  ON  ERROR  GOTO  2550:  GOTO  2560 
2550  RESUME  2530 
2560  OPEN  "O",  #1,  D$ 

2570  WRITE  #1,  EAD,  DIP,  OOND,  ELM,  ETM,  EVM,  DT,  IT,  AIT 

2580  WRITE  #1,  NS,  CE,  VEK,  AE,  CT,  VTK,  AT,  Z,  C,  WK,  EM,  CMLN,  CM,  FAR,  PF, 
SIG,  ST 

2590  WRITE  #1,  IRM,  DS,  NC,  E$,  P$,  K$,  C$ 

2600  FOR  1  =  0  TO  NC  —  1 

2610  WRITE  #1,  XO(I),  PDCC(I) ,  PDSL(I) ,  K(I) ,  E(I) ,  RO(I) 

2620  NEXT  I 

2630  FOR  I  =  0  TO  NC  -  1 
2640  FUR  J  =  0  TO  NS  -  1 
2650  WRITE  #1,  ES(I,  J) 

2660  NEXT  J 
2670  NEXT  I 
2680  CLOSE 
2690  GOTO  2860 

2700  INPUT  "program  data  file  name";  D$ 

2710  ON  ERROR  GOTO  2720:  GOTO  2730 
2720  RESUME  2700 
2730  OPEN  "I",  #1,  D$ 

2740  INPUT  #1,  EAD,  DIP,  OOND,  ELM,  ELM,  EVM,  DT,  IT,  AIT 

2750  INPUT  #1,  NS,  CE,  VEK,  AE,  CT,  VTK,  AT,  Z,  C,  WK,  EM,  CMLN,  CM,  FAR,  PF, 
SIG,  ST 

2760  INPUT  #1,  LRM,  DS,  NC,  E$,  P$,  K$,  C$ 

2770  FOR  I  =  0  TO  NC  -  1 

2780  INPUT  #1,  XO(I),  PDCC(I) ,  PDSL(I) ,  K(I) ,  E(I),  RO(I) 

2790  NEXT  I 

2800  FUR  I  =  0  TO  NC  -  1 
2810  FOR  J  =  0  TO  NS  -  1 
2820  INPUT  #1,  ES(I,  J) 

2830  NEXT  J 
2840  NEXT  I 
2850  CLOSE 

2860  PRINT  :  A$  =  "a" 

2870  INPUT  "print  encounter  parameter  values  (y/n)";  A$ 
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2880 

2890 

2900 

2910 

2920 

2930 

2940 

2950 

2960 

2970 

EAD 

2980 

DIP 

2990 

OCND 

3000 

ELM 

3010 

EIM 

3020 

EVM 

3030 

DT 

3040 

IT 

3050 

AIT 

3060 

NS 

3070 

CE 

3080 

VEK 

3090 

AE 

3100 

CT 

3110 

VTK 

3120 

AT 

3130 

C 

3140 

WK 

3150 

Z 

3160 

EM 

3170 


IF  A$  =  "Y"  OR  A$  =  "y"  THEN  GOTO  2900 

IF  A$  =  "N"  OR  A$  =  "n"  THEN  GOTO  3290  ELSE  GOTO  2870 


IPRINT 

IPKENT 

IPKENT 

IPRINT 

IPKENT 

IPKENT 

LEKENT 

LPRINT 

LPRINT 

IPKENT 

IPRINT 

LPRINT 

LEKENT 

IPRINT 

IPRINT 

IPRINT 

IPRINT 

IPRINT 

IPRINT 

IPRINT 

IPRINT 

IPRINT 

IPRINT 

IPRINT 

IPRINT 

LPRINT 

IPRINT 

IPRINT 


"program  file  name 
"program  data  file  name 
"electric  data  file  name 
'•processing  data  file  name 
"kinematic  data  file  name 

"combined  electric,  processing  &  kinematic  data  file  name 
"electrometer  axis  direction  (decimal  degrees)  " 

"electrometer  axis  depression  angle  (decimal  degrees)  " 

"encounter  medium  conductivity  (mhos/meter)  " 

"electric  longitudinal  moment  (ampere-meters)  " 

"electric  transverse  moment  (ampere-meters)  " 

"electric  vertical  moment  (ampere-meters)  " 

"sampling  period  (seconds)  " 

"integration  time  (seconds)  " 

"adjusted  integration  time  (seconds)  " 


"number  of  samples  per  encounter 
"electrometer  course  (decimal  degrees) 
"electrometer  speed  (knots) 

"electrometer  altitude  (meters) 

"target  course  (decimal  degrees) 

"target  speed  (knots) 

"target  depth  (meters) 

"electrometer  relative  course  (decimal  degrees) 
"electrometer  relative  speed  (knots) 
"electrometer-target  vertical  separation  (meters) 
"electric  dipole  moment  (ampere-meters) 

"dipole  moment  azimuth  (decimal  degrees) 


M 


'  +  N$ 

'  +  D$ 

'  +  E$ 

'  +  P$ 

'  +  K$ 

'  +  C$ 
SPC(2) 

SPC(2) 

SPC(2) 

SPC(2) 

SPC(2) 

SPC(2) 

SPC(2) 

SPC(2) 

SPC(2) 

SPC(2) 

SPC(2) 

SFC(2) 

SPC(2) 

SPC(2) 

SPC(2) 

SFC(2) 

SPC(2) 

SPC(2) 

SPC(2) 

SPC(2) 

SPC(2) 


ti 


3180  IFRINT  "dipole  moment  depression  angle  (decimal  degrees)  SPC(2) ; 

CM 

3190  IFRINT  "distance  between  samples  on  the  relative  track  (meters)";  SPC(2) ; 
D6 

3200  IFRINT  "false  alarm  rate  (false  alarms  per  hour)  SPC(2)  ? 

FAR 

3210  IFRINT  "false  alarm  probability  ";  SPC(2)  ; 

PF 

3220  LfRINT  "electric  noise  (volts/meter)  ";  SPC(2); 

SIG 

3230  LFRINT  "maximum  lateral  range  (meters)  SPC(2) ; 

UW 

3240  IFRINT  "lateral  range  step  (meters)  SPC(2); 

ST 

3250  IFRINT  "number  of  lateral  range  function  values  ";  SPC(2)  ; 


NC 

3260  FOR  I  =  0  TO  30 
3270  IFRINT 
3280  NEXT  I 
3290  PRINT  :  A$  =  "a" 

3300  INPUT  "print  lateral  range  function  values  (y/n)";  A$ 

3310  IF  AS  =  "Y"  OR  A$  =  "y"  THEN  GOTO  3330 

3320  IF  A$  =  "N"  CR  A$  =  "n"  THEN  GOTO  3410  ELSE  GOTO  3300 

3330  IFRINT  D$;  "  lateral  range  function  values" 

3340  IFRINT  :  IFRINT 

3350  IFRINT  "L  p(oc)  p(sl)  E  R 

k" 

3360  IFRINT  "meters  volts/meter  meters" 

3370  IFRINT 

3380  P0RI  =  0TDNC-1 

3390  IFRINT  X0(I);  TAB(IO) ;  PDOC(I) ;  TAB(24) ;  PDSL(I) ;  TAB(38)  ;  E(I) ;  TAB(52) ; 
R0(I);  TAB (70);  K(I) 

3400  NEXT  I 

3410  IRINT  :  A$  =  "a" 

3420  INPUT  "print  electric  signal  values  (y/n)";  A$ 

3430  IF  A$  =  "Y"  OR  AS  =  "y"  THEN  GOTO  3450 

3440  IF  A$  =  "N"  OR  AS  =  "n"  THEN  GOTO  3630  ELSE  GOTO  3420 

3450  PRINT  ;  PRINT  "identify  signed  by  encounter  lateral  range  index" 

3460  PRINT  "lateral  range  equals  the  index  times  "  +  STR$(ST)  +  "  meters" 

3470  PRINT  "index  values:  +  STR$((NC  -  1)  /  2)  +  "  to  "  +  STR$((NC  -  1)  / 

2) 

3480  INPUT  "lateral  range  index";  K;  I  =  K  +  (NC  -  1)  /  2 

3490  IFRINT  :  IFRINT  :  IFRINT 

3500  IFRINT  D$;  "  electric  signed  values" 

3510  IFRINT  "for  a  lateral  range  of  ";  K  *  ST;  "  meters” 

3520  IFRINT 

3530  IFRINT  "relative  CPA  distenoe";  TAB (35) ;  "electric  signed" 

3540  IFRINT  "s  in  meters";  TAB (35) ;  "Es  in  volts/meter" 

3550  IFRINT 

3560  FCR  J  =  0  TO  NS  -  1 

3570  IFRINT  (J  -  (NS  -  1)  /  2)  *  DS;  TAB (35) ;  ES(I,  J) 

3580  NEXT  J 


3590  PRINT  :  A$  =  "a” 

3600  INPUT  "print  electric  signal  values  for  a  different  lateral  range  (y/n) ' 
A$ 

3610  IF  A$  =  "Y"  CR  A$  =  "y"  THEN  GOTO  3450 

3620  IF  A$  =  "N”  OR  A$  =  "n"  THEN  GOTO  3630  ELSE  GOTO  3600 

3630  PRINT  :  A$  =  "a" 

3640  INPUT  "continue  to  vise  the  program  (y/n)";  A$ 

3650  IF  A$  =  "Y"  CR  A$  =  "y"  1HEN  Gai°  10 

3660  IF  A$  -  "N"  CR  A$  =  "n"  THEN  GOTO  3670  ELSE  GOTO  3640 

3670  END 

3680  R  =  SQR(X  *  X  +  Y  *  Y) :  REM  rectangular  to  polar  conversion 
3690  IF  R  =  0  THEN  T  =  0:  RETURN 

3700  IF  ABS (X  /  R)  =  1  THEN  Q  =  SGN(X)  *  (FT  /  2)  ELSE  Q  =  ATN(X  /  R  /  SQR(1 
X  *  X  /  R  /  R) ) 

3710  IF  ABS(Y  /  R)  =  1  THEN  T  =  (PI  /  2)  *  (1  -  SGN(Y) )  EISE  T  =  (PI  /  2)  - 
MN(Y  /  R  /  SQR(1  -Y*Y/R/R)) 

3720  IF  Q  <  0  THEN  T  =  2  *  PI  -  T 
3730  RETURN 

3740  Y1  -  XI:  IF  XI  <  0  THEN  Y1  =  -Yl:  REM  normal  approximation 
3750  G  =  1  /  (1  +  QO  *  Yl) 

3760  Yl  =  EXP(-Y1  *  Yl  /  2)  /  SQR(2  *  PI)  *  G  *  (Q1  +  G  *  (Q2  +  G  *  (Q3  +  G  ' 
(Q4  +  G  *  Q5) ) ) ) 

3770  IF  XI  >=  0  THEN  Yl  =  1  -  Yl 
3780  RETURN 
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